Direct Measurement of the Zak phase in Topological Bloch Bands 
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Geometric phases that characterize the topological properties of Bloch bands play a fundamental 
role in the modern band theory of solids. Here we report on the direct measurement of the geometric 
phase acquired by cold atoms moving in one-dimensional optical lattices. Using a combination of 
Bloch oscillations and Ramsey interferometry, we extract the Zak phase - the Berry phase acquired 
during an adiabatic motion of a particle across the Brillouin zone - which can be viewed as an 
invariant characterizing the topological properties of the band. For a dimerized optical lattice, 
which models polyacetylene, we measure a difference of the Zak phase equal to <5</?zak = 0.97(2)7r 
for the two possible polyacetylene phases with different dimerization. This indicates that the two 
dimerized phases belong to different topological classes, such that for a filled band, domain walls 
have fractional quantum numbers. Our work establishes a new general approach for probing the 
topological structure of Bloch bands in optical lattices. 



The non-trivial topological structure of Bloch bands 
in solids gives rise to fundamental physical phenomena, 
including fermion number fractionalization [IH3], the 
quantum Hall effect [U |5J, as well as topologically 
protected surface states in topological insulators [SJ \7\. 
The topological character of a Bloch band is defined 
by certain invariants, which can be expressed in terms 
of the Berry's phase [5] acquired by a particle during 
adiabatic motion through the band [SJ |7j- The most 
well-known example is the two-dimensional topological 
invariant, the first Chern number, which is related to 
the Berry's phase for a contour enclosing the Brillouin 
zone and determines the quantized value of the Hall 
conductivity of a filled two-dimensional band [H [S]. 
For one-dimensional systems, topological invariants of 
Bloch bands have been discussed theoretically 0011], 
however never been measured in any experiment. 

Here we present direct measurements of Berry's phase 
and topological invariants of one-dimensional periodic 
potentials using systems of ultra-cold atoms in opti- 
cal lattices. Topological properties of one dimensional 
solids are characterized by the so-called Zak phase - 
the Berry's phase picked up by a particle moving across 
the Brillouin zone [5]. For a given Bloch wave il>k(%) 
with quasimomentum fc, the Zak phase can be conve- 
niently expressed through the cell-periodic Bloch func- 
tion Uk{x) — e~ lkx ipk(%)'- 

,G/2 

<y5zak = i / (uk\dk\u k ) dk, (1) 

J-G/2 

where G = 2-Kjd is the reciprocal lattice vector and d is 
the lattice period [§]■ Non-trivial Zak phases underlie 
the existence of protected edge states [TOj HI] , fermion 
number fractionalization [IH3], and irrationally charged 
domain walls [HI Q3] between topologically distinct one- 



dimensional solids. These phenomena, initially discussed 
in the context of quantum field theory [TJ [21 [12] , later on 
found condensed matter realizations in polyacetylene [3J , 
described by the celebrated Su-Schrieffer-Heeger (SSH) 
model, and linearly conjugated diatomic polymers [13j . 

In our experiment, the key idea is to combine co- 
herent Bloch oscillations with Ramsey interferometry 
to determine the geometrical Zak phase and reveal the 
underlying topological character of the Bloch bands. 
Previously, the measurement of topological invariants 
was confined to two-dimensional bands by exploiting 
the relation between the Chern number and the Hall 
conductivity for a filled band introduced by Thouless- 
Kohmoto-Nightingale-DeNijs [4 . In contrast, in our 
case the integration over the Brillouin zone necessary for 
extracting topological invariants is achieved by adiabatic 
transport of a single-particle wave packet through the 
band using Bloch oscillations. Recently, it has also 
been suggested that in the context of ultra-cold atoms, 
topological properties could be studied through time- 
of-flight images [T4HT6] or measurements of anomalous 
velocity [5]IT7]. 

In the following we focus on a dimerized optical lattice 
with two sites per unit cell - a system which, despite its 
simplicity, exhibits rich topological physics, and depend- 
ing on the parameter values can mimic either polyacety- 
lene [3J, or conjugated diatomic polymers [T3]. Within 
a tight-binding model, the physics of such a system is 
captured by the Rice-Mele Hamiltonian |13) : 
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FIG. 1: Energy bands and topology of dimerized lattice 
model, (a) Schematic illustration of optical superlattice po- 
tential used in the experiment to realize the Su-Schrieffer- 
Heeger model (yellow box denotes the unit cell of size d = A s ). 

(b) Exemplary curves for the lower and upper energy bands 
(red and blue lines) and phase 6 k for dimerization Dl and D2 
(solid and dashed line) as a function of quasimomentum k. 

(c) Pseudo-spin representation of the eigenstates u T ,fe of the 
upper and lower energy bands for the two dimerization con- 
figurations Dl and D2. The pseudo-spin vectors u^k point 
in opposite directions and exhibit the same sense of rotation 
(winding) with quasimomentum k. In the phase Dl (D2) 
u T ,k evolve (anti-) clockwise and therefore exhibit opposite 
winding. 



Dl and D2 features fractionalized excitations. When the 
on-site energies are tuned to be different (A ^ 0), our 
system models a linearly conjugated diatomic polymer; 
in this case, the difference of the Zak phases is fractional 
in units of 7r, which generally gives rise to irrationally 
charged domain walls [UJ ■ 

In the experiment, we realized the Hamiltonian H of 
Eq. |2| by loading a Bose-Einstein condensate of 87 Rb 
into a one-dimensional optical superlattice potential |18) . 
This potential was formed by superimposing two stand- 
ing optical waves of wavelengths A s = 767 nm and 
A; = 2A S = 1534 nm that generate a lattice potential of 
the form V{x) = Visin 2 (fc;x + cj>/2) + V,sin 2 (2fc ; a; + tt/2), 
where fc; = 2it/Xi (Fig. [l^,). Phase control between 
the two standing wave fields enabled us to fully control 
<f>. For example, switching between = and <p = tt 
allowed us to rapidly access the two different dimerized 
configurations Dl (D2) with A = in the experiment, 
whereas by tuning 4> slightly away from these symmetry 
points, we could introduce a controlled energy offset A. 

The eigenstates of H can be written as Bloch waves of 
the form: 

ip k (x) = e lkx u k (x) = ^2 a k e lkXn w a {x - x n ) 

n 

+ W (x " +<i/2 W - £n - d/2), 

where k denotes the quasimomentum, x n = nd with n 
integer, and w a j,(x) are the Wannier functions [5U] for 
a n ,b n sites, respectively. The coefficients a k ,/3 k are de- 
termined through the eigenvalue equation Htpu = E k ip k . 
In this case, the cell-periodic wave function u k can be 
viewed as a two-component spinor = (a k , /3k), and 
Eq. ([I]) for the Zak phase takes an especially simple form: 



VZak = I 
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where J, J' denote modulated tunneling amplitudes 
within the unit cell, 6^(6^) are the particle creation 
operators for an atom on the sublattice site a n (b n ) in 
the nth lattice cell (Fig. [iji), and A characterizes the 
energy offset between neighboring lattice sites. 

When the on-site energies of the two sites are tuned 
to be equal (A = 0), our system corresponds to the 
SSH model of polyacetylene (see Fig. [iji). In this 
case, H is known to exhibit two topologically distinct 
phases, Dl for J > J', and D2 for J < J', separated 
by a topological phase transition point at J = J'. 
The distinct topological character of the two phases is 
reflected in the difference of their Zak phases, for which 
<5</?Zak = 7T- At half filling, a domain wall between phases 



For our choice of the unit cell [35], the eigenfunctions 
for the lower (upper) band of the SSH model (A = 0) are 
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where 9 k is determined through J e lkd / 2 -\- j' e -ikd/2 _ 
|£fc|e* efc (see Supplementary Information). We can thus 
visualize the Bloch periodic functions as pseudo spin- 1/2 
states oriented in the equatorial plane of a Bloch sphere 
(Fig. [I];). Note also that although ip k+ G(x) — tpk(x), this 
translational invariance is not true for it^fe, because in 
our system with a two-site unit cell u T . k +a = ^u^^, 
where a z is the third Pauli matrix. As the two state 
vectors for the upper and lower bands are orthogonal, 
they point in opposite directions and therefore exhibit 
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FIG. 2: Experimental sequence and spin-dependent Bloch os- 
cillations, (a) Energy band, MW pulses and state evolution 
of a single atom in a superposition of two spin-states with op- 
posite magnetic moment (brown and green balls) during the 
three-step echo sequence described in the text. The winding 
of the state vector with k is given by Ok (solid line dimeriza- 
tion Dl, dashed line dimerization D2). (b,c) Time-of-flight 
momentum distributions taken for different evolution times of 
the spin-dependent Bloch oscillations in the lower (b) and up- 
per energy band (c) used in the experiment. Each momentum 
point is an average of three identical measurements. 



the same winding when the quasimomentum k is var- 
ied adiabatically. The Zak phases for the lower and up- 
per band are thus identical y^ak = ""/^ However, when 
the dimerization is changed from configuration Dl to D2 
(Fig. , the corresponding geometric phase changes to 
fzak ~ — because of the opposite winding of the 
state with quasimomentum k. The difference of the two 
Zak phases for the two dimerized configurations is then: 



<Vzak = ¥>Zak - <Pz 



ak 



(3) 



We point out that the Zak phase of each dimerization is 
a gauge dependent quantity, i.e. it depends on the choice 
of origin of the unit cell, however, the difference of Zak 
phases of the two dimerizations is uniquely defined [71l21j. 

When an atom is adiabatically evolved through the 
Brillouin zone of the periodic potential k — > k + G, it 
acquires a phase shift due to three distinct contribu- 
tions: i) a geometric phase tpzak as well as ii) a dynamical 
phase <pdyn = f E(t)/hdt, both derived from the band- 
structure, and iii) a phase due to the Zeeman energy of 



the atom in an external magnetic field (see semiclassical 
analysis in Supplementary Information): 

<Ptot — <PZak + fdyn + fZeeman- 

To isolate the geometrical Zak phase in the exper- 
iment, we employ a three-step sequence (Fig. [2^, and 
Supplementary Information). Step 1) We start with an 
atom in the state k = 0) and bring it into a coherent 
superposition state 1/V2(|t; k = 0) + ||, k = 0)) using 
a microwave 7r/2-pulsc. Here a =t,| denote two spin 
states of the atom with opposite magnetic moment. 
Then a magnetic field gradient is applied that creates 
a constant force in opposite directions for the two spin 
components. Such a constant force leads to Bloch oscil- 
lations, i.e. a linear evolution of quasimomentum over 
time [TH]- In our case the force is directed in opposite 
directions for the two spin components. The atomic 
wavepacket thus evolves into the coherent superposition 
state l/V^Gt, fc) + e i5 v\\,-h)). When both reach the 
band edge, the differential phase between the two states 
is given by Sip — ipzak + <5^Zccman- Note that the dy- 
namical phase acquired during the adiabatic evolution is 
equal for the two spin states and therefore cancels in the 
phase difference. Step 2) To eliminate the Zeeman phase 
difference, we apply a spin-echo 7r-pulse at this point 
and also switch dimerization from Dl— >D2. For atoms 
located at the band edge k — ±G/2, this non-adiabatic 
dimerization switch induces a transition to the excited 
band of the SSH model. Step 3) The sequence is finally 
completed by letting the spin components further evolve 
in the upper band until they return to k = 0. At this 
point in time, a final 7r/2-pulse with phase <pmw is 
applied in order to interfere the two spin components and 
read out their relative phase Sip through the resulting 
Ramsey fringe. The change in dimerization occurring 
at the mid-point of the echo sequence is crucial in order 
not to cancel the Zak phase in addition to the Zeeman 
phase. Due to the opposite windings of the Bloch states 
in the upper and lower bands with quasimomentum k 
(Fig. [lp) , the resulting phase shift encoded in the Ramsey 
fringe is thus given by: Sip = V?zak — Vzak ^ the dimer- 
ization is swapped, whereas Sip = if it is left unchanged. 



In Fig. [2]d,c we show images of the momentum 
distribution of the atoms during the spin-dependent 
Bloch oscillations in the lower and upper energy bands. 
Note the opposite evolution in momentum space due to 
the opposite magnetic moments of the two spin-states. 
Atoms in the upper energy band are characterized by a 
distinctively different momentum pattern from atoms in 
the lower energy band. The Bloch oscillations period of 
iBioch = 0.85(3) ms was chosen to be slow enough, such 
that non-adiabatic Landau-Zener transitions at the band 
edge are negligible, while still maintaining an overall 
fast evolution time to minimize decoherence effects. 
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FIG. 3: Determination of the Zak phase, (a) Following the 
sequence described in the text, the atom number in the two 
spin states Nf 4 is measured and the fraction of atoms in the 
|t) spin state rif = Nf/(Nf + iVj.) is plotted as a function of 
the phase of the final microwave 7r/2-pulse. The difference in 
phase of the two Ramsey fringes yields the Zak phase differ- 
ence (Vzak = ip§*k — ipz&k- Blue (black) circles correspond 
to the fringe in which the dimerization was (not) swapped. 
In order to reduce the effect of fluctuations, every data point 
is an average of five individual measurements and the error 
bars show the standard deviation of the mean. The phase of 
the reference fringe (black) is determined by a small detuning 
of the microwave pulse [32] . (b) Measured relative phase for 
14 identical experimental runs (left), which give an average 
value of Sipzak = 0.97(2)-7r. The corresponding histogram is 
shown on the right with a binning of 0.05-7T. The lcr-width of 
the resulting distribution is a = 0.077T. 



A typical result for the two Ramsey fringes obtained 
with and without dimerization swapping during the state 
evolution can be seen in Fig. [3^,. Each plotted value for 
a given angle <pmw is an average over five identical mea- 
surements in order to reduce the effect of residual fluc- 
tuations. We performed a further statistical analysis by 
recording 14 independent Ramsey fringes for the two con- 
figurations. The obtained phase differences are shown in 
Fig. [3]d together with the corresponding histogram. From 
these individual measurements we determine the geomet- 
ric phase difference between the two dimerized configu- 
rations to be: 




Sep = 0.97(2)71-, 



FIG. 4: Fractional Zak phase, (a) Lattice potential with- 
out and with an on-site energy staggering A. When A = 
the Zak phase is <pzak(A = 0) = tt/2. As A increases, the 
pseudo-spin vectors move away from the equatorial plane and 
the value of y>zak(A) decays rapidly to zero, (b) Measured 
phase difference (^Zak — (^zak(A) as a function of A. Each 
individual point was obtained from four individual measure- 
ments. The vertical error bars represent the standard error of 
the mean. The green line is the theoretical prediction and the 
shaded area represents the uncertainties in the calibration of 
the energy offset A. The insets show a typical Ramsey fringe 
for A/ J = -1.2 (left) and A/J = 1.2 (right), which were used 
to extract the relative phase Sep. The blue (black) fringes cor- 
respond to measurements with (without) staggering |32| . 



in excellent agreement with theory, as discussed above. 
The uncertainty in the recorded value denotes the 
standard error of the mean obtained from the distri- 
bution function (Fig. [3)3) and is mainly determined 
by experimental imperfections in the control of the 
underlying lattice potentials. 

To further investigate the variation of the Zak phase 
with lattice parameters, we performed a second series 
of experiments with a staggered on-site energy offset 
A (Fig. [i^), corresponding to a heteropolar dimer 
configuration [13]. The energy offset A displaces the 
pseudo-spin Bloch vectors away from the equatorial 
plane, resulting in an additional dependence of the 
Zak phase on the offset A (Fig. [4^, and Supplementary 
Information). In order to probe the dependence of ty?Zak 
on A, we performed an experimental sequence that was 
similar to the one described above. However, instead of 
swapping the dimerization from Dl to D2, an energy 
offset |A| < 2 J was introduced for one half of the 
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sequence. Thereafter, because of the spin-echo pulse, 
the wavepackets return to k — in the lowest band. 
Though the system completes a full Bloch oscillation in 
the lowest band, the total geometric phase acquired is 
not zero, since the Bloch vector is displaced from the 
equatorial plane during one half of the sequence (Fig. |4| 
and the Zak phase is changed from ipzak to <,5zak(A). 
The resulting phase in the Ramsey fringe is thus given by 
Sip = ipz&k ~~ <^Zak(A) when the energy offset is present, 
and it is Sip = when the offset is absent. As before, the 
phase difference between these two fringes for atoms in 
the lowest band allows us to determine the relative phase 
</?Zak — ^Zak(A). During the non-adiabatic switching of 
the superlattice potential at step 2 of the experimental 
sequence some of the atoms are transferred to the higher 
band and acquire a different geometric phase. However, 
taking into account this contribution to the measured 
phase difference enabled us to extract the relative phase 
Sip from our data (see Supplementary Information). As 
shown in Fig. ^p, we find good agreement between the 
measured and predicted values of the fractional Zak 
phase. 

In conclusion, we have presented a general approach 
for studying topological properties of Bloch bands in 
optical lattices and demonstrated its versatility through 
a first direct measurement of the topological invariant 
in topologically non-trivial Bloch bands. Making use of 
the recently demonstrated control of optical potentials 
at the single-site level [22], we plan to realize domain 
walls in the dimerizcd lattice that would allow us to 
directly study edge states [33] IMj and fractional charges 
located at the interface of the two topologically distinct 
phases [TH21 H5J. Although in this work we focused 
on one-dimensional systems, our technique can easily 
be extended to two-dimensional systems, where the 
change of the Zak phase in the Brillouin zone gives the 
topological density of the Bloch band [5S] . This enables 
measurements of both the Chern number of topological 
bands and the 7r-flux associated with a Dirac point. 
Additionally, we expect that this idea can be extended to 
measure the non-Abelian Berry's phase in Bloch bands, 
such as in a system with quantum spin Hall effect |27j . 
to the study of Floquet states in periodically driven 
systems [28-30], and to quasiparticles in unconventional 
superconductors, such as d-wave superconductors, which 
have Dirac dispersion at the nodal points [31]. Overall, 
our work indicates that cold atomic systems provide 
a versatile platform for studying topological states of 
matter, and establishes a novel method for probing their 
properties. 
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In this section we discuss the Bloch wave functions 
in the SSH model, and calculate the Zak phase of two 
dimerized configurations of the model. We consider the 
most general case when not only tunneling but also on- 
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site potentials on sublattices A and B are dimerized |S1| : 
H = -J2{ M nk + J'ajA-i+h.c.) (SI) 

n 

+ A X^Ja- t>l k)- 

71 

The energy bands and eigenstates can be obtained from 
the eigenvalue equation 



A -p fc 
-Pi -A 



a k \ ~ a k 



where 



p fc = Je lkd/2 + J' e - lkd / 2 = \e k \e l9 K 



(S2) 



(S3) 



The variable e k in the above expression corresponds to 
the energy of Bloch states in the SSH model (A = 0). 



Eigenstates 



Solving Eq. (S2|, we obtain for the eigenenergies 



(S4) 



Following Ref. [S2| . we require that a k +c = a k and 
fik+G — —flk- The eigenstates in the lower and upper 
bands which satisfy this condition can be chosen as fol- 
lows: 



(S5) 



where 7fc and k are given by (the expression for 9 k is 
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^ sin^e"^ 



obtained from Eq. < S3 1 ) 



7fc = arctan 



A 



arctan 



(J — J') sin(fcd/2) 
(J + J')cos(kd/2)' 

Zak phase 

Generally, the Zak phase for a given band is given by 

r G/2 



<PZak = I I (a* k d k a k + /3* k d k /3 k )dk. 

I-G/2 



For the wave functions in Eq. (S5|, this gives 

r G/2 



</5 Zak - (A) 



-G/2 
G/2 



cos 2 y d k 9 k dk 7 



<p Z ak,+ (A) = / sin 2 ^-d k 9 k dk 

G/2 1 



(S6) 



(S7) 



These equations were used to make the green theoretical 
curve in Fig. 4 of the main text. 



For the SSH model A 

2 , 



such that cos 2 7fc/2 = 
sin 2 7fc/2 = 1/2 and the Zak phases reduce to fzlk - = f 



ZAK PHASE AND THE CHOICE OF A UNIT 
CELL 

We emphasize that the Zak phase, in general, depends 
on the choice of unit cell for the crystal. Here we would 
like to clarify why for our choice of the unit cell of the 
SSH model the Zak phases are given by ±7r/2 (rather 
than by or 7r, as in some previous papers, e.g., see 
Ref. jS3| and references therein). 

To understand the relation of the Zak phase to the 
choice of the unit cell, we first note that the function 
u k (x) is not periodic under a translation by a reciprocal 
lattice vector G = 2it/d, but rather satisfies the following 
condition (see Ref. [52]): 



u k+G (x) = e lGx u k (x) 



(S8) 



This guarantees that the full Bloch function 
ipk(x) = e lkx u k (x) is periodic under a translation 
by a reciprocal lattice vector G. 

The Zak phase is defined in terms of the periodic Bloch 
functions (see Ref. |S2j): 



"PZak 



G/2 



G/2 



A(k)dk, A(k) = i(u k \d k u k ) 



(S9) 



where A{k) is the Berry's connection, and the scalar 
product is defined as (u\v) = u* {x)v{x)dx. 

Let us now consider a translation by a distance a: x' = 
x + a, which leads to a different choice of the unit cell 
(unless a = nd, where n is integer). The new unit cell 
is x 1 £ [0,d), or, cquivalently, x G [— a, d — a). Then, 
in order for the condition in Eq. ( S8 1 to be satisfied, the 



Bloch functions have to be redefined as follows: 

u' k {x') = u k {x' - d)e- lka . (S10) 

The Berry connection is modified accordingly: 

A'{k) = {u' k \d k u' k ) = (u k \d k u k ) -ia= A(k) - ia. 

Integrating the Berry connection over the Brillouin 
zone, we obtain the change in the Zak phase: 



'y'zak = ^Zak + Ga = ipz + 2ira/d. 



(Sll) 



The Zak phase changes by 2tt under a translation by a 
multiple of the lattice vector d. Thus it is not surprising 
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that our choice of the unit cell for the SSH model gives 
rise to the values of the Zak phase which are different 
from the values mentioned in some other works [S3] . 
The difference of tt/2 appears due to the shift of the unit 
cell by a = d/i. 

We also emphasize that, even though the Zak phase 
itself depends on the choice of the unit cell, the differ- 
ence of the Zak phases in two states, measured in our 
experiment, is an invariant which describes the topolog- 
ical properties of the two states (see Ref. |S4|). 



BLOCH OSCILLATIONS IN SUPERLATTICES 



and 



[H,, 



,q c — ,q 



F 
N 
d 



^2 \a-,q(xn - x a )e iqXri al 



we find that the Ansatz in Eq. ( S13 1 provides a solution 



of the Heisenberg equation of motion when v = f = F/h 
and 



-iA =e_ ifco _ /t A + fx Q A 
fB 



fA 



(u-,fc n -/t|<9fc M -,fco-/t) 



(u- jko -ft\d k u +lko -f t ) 



(S15) 



Here we derive the dynamical equations for the Bloch 
oscillations in the dimerized lattice. We consider a dimer- 
ized lattice subject to an external force, described by the 
Hamiltonian 

H F = H - F | (x„ - ar )a^d n + (x n + ^ - x )P n b n | 

with x n = nd. For simplicity, we wrote this equation 
for only one of the pseudo-spin species; an equation 
for the dynamics of other pseudo-spin is obtained by 
a substitution F — > — F in the above equation. Slowly 
fluctuating magnetic fields in our experiment will add 
an extra Zeeman energy difference between two pseudo- 
spins, which can be absorbed in the definition of xq. 

Eigenstates of Hp are plane waves with energies ztik — 
±-\/A 2 + e\. For our analysis it is convenient to work in 
second quantization, 

v n 



and a similar equation for B. The last term in Eq. (S15) 



e lkx "al 



(S12) 



where N is the total number of sites in the lattice. For 
finite F we need to solve the Heisenberg equation of mo- 
tion f t &{t) = \[H F ,&(t)]. We look for solutions where 
the quasimomentum changes at a constant rate, 



t 

-\-,ko — vt' 



(S13) 



describes non-adiabatic mixing of the bands, which we 
neglect. Assuming that atoms occupy only the lower 
band (B = and \A\ = 1) and taking A(t) = e**®, 
we obtain: 

/ 

^ = e-,k -ft + fxo - j{u- ko -f t \d k u- t k -ft). (S16) 



The first term in Eq. (Sf6) describes the dynami 



cal phase contribution, the second is, effectively, the 
Zeeman phase (which contains the effect of the fluc- 
tuating magnetic fields that shift xo), and the last 
term describes the Berry's phase part. Integrating 



Using 



Eq. (S16l over a period of the Bloch oscillations 
gives Eq. (3) of the main text, with yzeeman = ±fxot, 
where ± sign corresponds to the two pseudo-spin species. 

A crucial point is that our spin-echo type protocol is 
insensitive to both the Zeeman phase (or, equivalently, 
to the value of xo) and the dynamical phase. The spin- 
echo nature of the experiment guarantees that the Zee- 
man phases for the two trajectories are equal (each atom 
spends half of the time in the spin-up state, and half of 
the time in the spin-down state). The dynamical phases 
acquired during the motion in the Bloch bands cancel out 
due to the reflection symmetry of the dispersion relation, 
Ek = E-k- Both of these points finally allow us to single 
out the Zak phase. 



EXPERIMENTAL SEQUENCE 



dt 



E 



'N 

d/3-, k 
dk 



E 



10L— Jz&ri 

i(3-, k {x n 



da-, k 
dk 

4) 



e a' n \k=k a -vt 



e u n \k=k — vt 



(S14) 



The experimental sequence started by loading in 
200 ms a Bose-Einstcin condensate of about 5 x I0 4 
atoms in the Zeeman state \F — l,mp = —I) = |^) into 
a dimerized ID lattice of Vi = l3E r j, V s = 4E r , s and 
(f> = 0, where E r ^ — h 2 /2m\i, i — Z,s. Thereafter, a 
magnetic field gradient of / = 1.18(4) kHz/d was ramped 
up in 1.5 ms, and subsequently a 3-fj.s MW 7r/2-pulse 
coupling to the Zeeman state \F = 2,m F = —1) = |t) 
was applied. After half a Bloch oscillation, when the 
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First long lattice, t(> = 



Second long lattice, ■ 



Long Lattices 



V(x) 



WWVW VWWW 




Superlattice 
potential 



J J' J J' J J 
Dimerization D1 



J J' J J' J J' 
Dimerization D2 



FIG. SI: Left (right): Superlattice potential created by using 
the first (second) long lattice with a phase cj> = (4>' = 7r). 
The corresponding dimerization is Dl (D2). 



atoms were at the edge of the band, a 7r-pulse was 
applied and the dimerization was swapped (See section 
Dimerization swap) . After an additional time of 425 /is 
the atoms returned to k = and a final 7r/2-pulse with 
a phase ifMW wa s applied. At the end of the sequence 
we measured the population fraction of each of the 
two spin components by using time-of-flight imaging 
with a Stern-Gerlach gradient field applied during the 
expansion to separate the different Zccman states. 

For the 850-/is Bloch oscillation time and 
J I J' = 1 kHz/0. 06 kHz, a negligible fraction of 
atoms was transferred to the upper band because of non- 
adiabatic Landau-Zener transitions at the band edge. 
This underlines the adiabaticity of the Bloch oscillations. 

Due to a reproducible 80-kHz drift in the offset mag- 
netic field during the time of a full Bloch oscillation, the 
three MW pulses were not all on resonance. The first 
(last) 7r/2-pulse was detuned by —40 kHz (40 kHz) and 
the intermediate 7r-pulse was on resonance, while the cou- 
pling strength is about 80 kHz. The drift in the offset 
magnetic field adds a constant value to the Zeeman phase 
^Zeeman, which is independent of the change in dimer- 
ization/staggering, therefore it does not affect the phase 
difference Sip. 
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FIG. S2: Atom fraction in the lower band (n_ ) as a function of 
the energy offset in units of J. As A/ J increases the fraction 
transferred to the upper band (n+) also increases. In the limit 
A/J — > co, n_ — > 0.5 



MEASURING THE ZAK PHASE AT A / 0: 
PROJECTION ONTO LOWER AND UPPER 
BANDS 

Our protocol for measuring the Zak phase at A / 0, 
described in the main text, involves a sudden turn on/off 
of A when the atoms, initially in the lower band, reach 
the edges of the Brillouin zone at k = ±ir/d. In such 
a process, a fraction of atoms get excited to the upper- 
band. The upper- and lower-band populations n± after 
such a non-adiabatic turn on/off of the staggering A can 
be obtained by projecting the lower-band eigenstates in 
Eq. (IS5} with A = onto the states with A/0. They 



are given by 



n± 



1 =F sin 7 fc _ £ fc T £fc 
2 ~ 2s k 



(S17) 



In the above expression, all quantities are defined for the 
final non-zero value of A. 



DIMERIZATION SWAP 

In order to change the dimerization from Dl to D2 
we employed a second long lattice, whose amplitude and 
relative phase </>' can be independently controlled. The 
dimerization exchange was performed by quickly switch- 
ing off the first long lattice with phase <\> — and a 
ramp-up of the second long lattice with a phase 4>' = tt 
within 10 /is. As the dimerization swap time is much 
shorter than the Bloch oscillation time, the dimerization 
exchange can be considered instantaneous. The resulting 
lattice potentials of the different lattice configurations are 
shown in Fig. |S1| corresponding to the two dimerizations 
of the SSH model. 



EXTRACTING THE PHASE DIFFERENCE 

V?zak - <^zak(A) FROM THE MEASURED RAMSEY 
FRINGES 

In this section we explain how we extract the value of 
^Zak~ <^Zak(A) from the measured fringes for Fig. 4 of the 
main text. Due to a fraction of atoms transferred to the 
upper band when A/0, the phase difference between 
the two Ramsey fringes deviates from <pzak — <Pzak(A). 
The reason for this is that the fraction of atoms in 
the upper band acquires a different geometric phase 
9Szak,+ (A) = tt - <yS Zak -(A) (see Eq. ([S7])) than the 
fraction in the lower band. Therefore, when measuring 
the Ramsey fringe, both phases enter into play as 
explained below. 
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FIG. S3: Measured phase difference tfExp as a function of the 
energy offset in units of J. The green line shows the predicted 
value of cpzak — (pzak(A) when there is no population transfer 
to the upper band. 



ffLi 



iZ|0}_< 



-i^zak,-(A)/2 



-iVZak,+ (A)/2 1 i^zak/2 



II) 



:|o). 



i VZafc ._(A)/2 
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Finally, we apply a 7r/2-pulse with a MW phase that 
rotates the spins to |f) -> -L( e -w«'/ 2 ||) + e +VM «'/ 2 ||)) 

and ||) -^{e VMW l 2 \\) +e- VMW / 2 \i)). After this pulse, 
the total population in |f) is given by: 



n t = 2 i 1 + n ~ C0s((pzak - (^Zak,-(A) + ¥>Afw) 

+ n+ cos((^ Zak - <^zak : +(A) + <^mw)} j (S18) 



As described in the main text, for half of the measured 
points during the first part of the Bloch oscillation, the 
atoms evolve in the lower band of the dimerized system 
with A = (for the second half of the measured points 
the atoms start in the lower band with A / 0, which 
gives essentially the same result). When they reach the 
edge of the band, the wave function is 

= -^(|t,-G/2)_e^ k / 2 + ||,G/2)_e^ z - k / 2 ), 
v2 

where <pzak — ¥>Zak,-(A = 0). When the MW 7r-pulse 
is applied and the energy offset A is quickly introduced, 
the spins are flipped and the population in the bands are 



n± (see Eq. (S17] and Fig. |S2|) 
1 / 



\, - = = - G/2)_ G/2) + )e^-/ 2 ||) 

1 / 



V2 V 

J^nZ\G/2)_ + V^|G/2)+)e- i ^/ 2 |t). 



In the second half of the Bloch oscillation the particles 
return to k = and the fraction in the upper (lower) 
band picks up a phase <^zak,+ (A) (<pzak,-(A)), such that 
the resulting state vector is: 



and a similar expression for n±. The Ramsey fringe 
is obtained by measuring the population in each spin 
component as a function of the MW phase ifMW- 

For negligible n + the phase of the Ramsey fringes 
directly corresponds to </?Zak — <pzak,-(A). For any fi- 



nite value of n + , we can recast Eq. ( S18 1 in the form 
A + B cos(ip Ex p + <£mw) and solve for t^zak - </?zak ,-(A) 
using the relation (/? Z ak,+ (A) = tt — <y2zak,-(A) and the 
theoretical values of n±, where (fExp is the measured 
phase difference between the Ramsey fringes. The re- 
sulting values of </?zak — ^Zak,-(A) are displayed in Fig. 4 
of the main text, and the values of if Exp are shown in 
Fig. [S3] 
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